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Abstract 

In this paper we sum over the spherical modes appearing in the 
expression for the Casimir energy of a conducting sphere and of a 
dielectric ball (assuming the same speed of light inside and outside), 
before doing the frequency integration. We derive closed integral ex- 
pressions that allow the calculations to be done to all orders, without 
the use of regularization procedures. The technique of mode summa- 
tion using a contour integral is critically examined. 



1 Introduction 



The Casimir effect is a remarkable consequence of QFT, exhibiting the reality 
of the zero point energy of the vacuum. This energy can be made manifest by 
studying its dependence on various constraints imposed on the electromag- 
netic field. Casimir's original analysis for parallel plates yielded a force 
per unit area {h = c = 1) 



240a4 ' ' 

between two conducting plates separated by a distance a. The Casimir force 
for a conducting spherical shell was first calculated by Boyer p|. The result 
was indeed puzzling: the force turned out to be repulsive. This was surprising 
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in view of the connection between Casimir's force and Van-Der-Waals forces 
which are attractive. This phenomenon may be also related to a fundamen- 
tal difference between the spherical Casimir problem and that of the plates: 
expansion of the sphere changes the distances between the points on the shell 
itself, and thus it is not clear weather the Casimir energy (associated with 
interaction among different patches on the shell) and electromagnetic the self 
energy of the shell can be separated It is not surprising, however, that 
this in turn called for a large number of subsequent calculations by various 
techniques @, ^ |, |^, §. There has been some controversy concerning the 
regularization procedures used in the calculations, bearing, for example on 
the relevance of Casimir's effect to sonoluminescence |lT| |12|, [13|. Since 
the repulsive nature of the Casimir force for spherical boundary conditions 
is somewhat counterintuitive it seems of importance to justify it more rigor- 
ously. 

We define the Casimir energy of a conducting sphere as follows: Let 
^(a,R) eigenmodes of the electromagnetic field inside two concentric 

spherical shells of radii a and R {a << R), where a large sphere of radius 
R is introduced as an infrared cutoff in order to discretize the spectrum 
outside the smaller sphere as well as inside it. The zero point energy of this 
system, namely | Yin ^rt'^^ is a divergent sum. The Casimir energy is then 
defined as the zero point energy difference between our configuration and a 
configuration where the inner sphere is of radius b {a << b < R) fl^ pl 



Ec = lim lim lim - V (e-"""''''''^"'^) - e-""" '""cui^'^)) (2) 

b->oo R^oo 2 ^ ^ n I \ I 

n 



Here e"'^'^" is an exponential regulator. We show that if the subtraction 
procedure is done properly, no further regularization (apart from keeping the 
naive exponential regulator, required to make the sums well defined, and 
the subtraction) is actually needed. In the sequel , we follow closely the 
procedure of mode summation using a contour integral, but our main result 
can also be applied to the Green's function method 

Most of the calculations have been carried out for the case of a conducting 
sphere and of a dielectric ball having the same speed of light inside and 



outside, that is, ^Jejl = v^e'/i' = 1 [|T5|, where e, are the permittivity and 



permeability of the medium surrounding a ball of radius a, and e', ii' are the 
permittivity and permeability inside the ball. It has been shown that the 



^this form of subtraction has the advantage of having a natural one to one correspon- 
dence between the eigenmodes of the constrained system (inner radius a) to those of the 
"unconstrained" one (where b and R are taken to oo) 
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Casimir energy in the latter case is given by [T^ 

1=1 

where x = ||/| , z/ = / + | and 

H + jj. a 



(3) 



(4) 



S/(Xj = lX]i{lX) 




2 e 2 J^[%x) 



(5) 



ei(x) = ixh[^\ix) = i 




Tree iTTv 
2'' 



Hl{ix) 



(6) 



and in the end of the calculation the regulator 5 is taken to zero. In the 
case ^ = 1, (3) is equivalent to the expression for the energy of a conducting 
sphere in vacuum. 

In order to evaluate (3), extensive use has been made of uniform asymp- 
totic approximations of the modified Bessel functions si{x) and ei{x). Our 
evaluation refrains from the use of these asymptotics, which, although very 
convenient in actual calculations, somewhat obscures the nature of the ex- 
pression (3). By studying the Z-sum of this expression before making the 
frequency integration, using simple orthogonality considerations we arrive at 
a general formula for sums of the form ^(2/-|-l)(s/e;)". We show that to first 
order in the integrand consists of a constant term (which is independent of 
the radius) and a function of a. Thus the divergence appearing in the eval- 
uation of (3) is due to integration of this constant over all frequencies. This 



divergence was already pointed out by Candelas [T^, however, this constant 
term would not be present if subtraction of the Casimir energy of a large 
radius sphere is done properly. 

We start by carefully reviewing the procedure of mode summation, using 
a contour integral, as was done in some previous works 0, |^, |T^. In this 
method the final integration is carried along the imaginary axis. Justification 
for the neglect of the part of the contour away from this axis will be given, by 
concrete estimates, taking into account the existence of an infinite number 
of poles of the integrand on the real axis. 
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2 Mode summation for conducting spherical 
shells 



In order to make the discussion well defined we investigate the following 
situation [^: A conducting sphere of radius a is placed (concentrically) 
inside a large conducting sphere of radius R. Using multipole expansion iST 



the appropriate modes of the electromagnetic fields can be found, through 
the boundary conditions 1?^ = and Eq^^ = on the shells. The eigenmodes 
are given as solutions to the equations: 
(I) Inside the inner sphere 



A 



(1) 



ji{uja) = 



(7) 



= -^[rji{ujr)] 



II) Between the shells 

Al = ji{ujR)ni{uja) - ji{uja)ni{uR) = 



(9) 



d 



d_ 

dr 



[rniiuor] 



r=R 



r=R 



d_ 

dr 



[rni{ujr)] 



= 
(10) 



Where equations (7), (9) refer to transverse electric modes, and (8), (10) refer 
to transverse magnetic modes. Next we define the function |jl9| 



Ai{uj;a,R) 



'11^ 



A; is an analytic function, whose positive real zeros are the modes of the 
electromagnetic field inside the concentric spheres. 

In order to sum over the modes of a given /, we use the identity [|^ Q 



N 

E 

n=l 



(n,l)i 



— (b dze~"''z-^\nAi{z 



2711 



dz 



(12) 



Here Cjy is a contour which encircles the zeros ^^(i ..ci;(jv/) of A^ . Since 
there is an infinite number of modes for each /, this sum will diverge if we 
take N to infinity and a to zero. As usual, the Casimir energy should be 
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Ln + ILn 




Figure 1: A rectangular contour in the complex plane. 



identified with that obtained after the subtraction of the zero point energy 



of the unconstrained "free" system ||14 



I °° I r 
AE(a, b, R) = lim V(/ + -) i dze^^^zFi 



(13) 



1=1 
Where 

Fi{z; a, 6, R) = In Ai{z; a, ^) - ^ In M^; R)) ■ (14) 
The Casimir energy is then given by the limit: 

Ec = lim lim lim AE{a,b,R). (15) 

After the limits R —>■ oo and b ^ oo are taken the remaining expression 
will have no reference to R and b. Removal of the ultra-violet cut off by 
taking the limit a — *• is done in the end of the derivation, after insuring 
convergence for all values of b and R. 

In order to make the evaluation of the contour integrals easier we take 
rectangular contours of side Ln (Fig. 1). Our first aim is to carry the entire 
integration along the imaginary axis, by showing that the contribution of 
the part of the contour away from it vanishes as L^v goes to infinity. The 
justification of this neglect, however, is not obvious, since the integrand in 
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(p!3| ) has infinitely many singularities along the real axis, which become ever 
denser as i? ^ oo. This follows from the property that for any finite R there 
is an infinite number of field modes, and subsequently there is a corresponding 
infinite number of poles of z-^lia Ai{z,b, R) along the real axis. Even after 
taking the limit R ^ oo the spectrum will still consist of a discrete set 
(eigenmodes inside the inner shell) imbedded in a continuum part. Indeed, 
let us examine the behavior of Ai{uj; a, R) for large arguments u. Using the 



asymptotic behavior pO 



1 In 

Ji - sm(x - y ) (16) 
1 , In 

ni cos(x - — ) (17) 

X 2 



we obtain 



Ai(uj; a, R) — — — sin(2ci;a — In) sin^(ci;(i? — a)) (18) 
za^Ru 

Thus for large arguments the modes behave as ^ or 2(R-a) where n is an 
integer. (Note that the density of modes for large frequencies is similar to 
that appearing in the one dimensional version of the standard Casimir effect 
for parallel conducting plates separated by a distance a inside a large box 
of length R.) Let us examine the behavior of the integrand for large values 
of \uj\. Substituting the asymptotics (|T^) , we show that the integral on a 
contour that does not pass through one of the poles, decays exponentially 
fast. The integrand is then simply: 

Fi{lj; a, b, R) = 2(a(cot(2co'a — In) — cot(a;(i? — a))) + 
R{cot{uj{R - a)) - cot{uj{R - b))) + 

b{cot{iu{R-a))-cot{2ujb-ln))) (19) 

It is easy to check that for Im — oo and any real A, B, 

I cot{Auj) - cot{Buj) I = o(e-2™"(^'^)^'^ '^) (20) 

Thus on the upper and lower parts of the contour (Im z = ±Lj\f) the integral 
is bounded by: 



AT 



{iL^ + z)Fi{iLN + z)dz 







JV 



< / \iLN + z\12Re-^''^^dz 
Jo 

< ISLlfRe"*"^'' (21) 
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The part of the contour with Re z = Ln can be evaluated as follows. First 
we prove that integrals of the form 



-a{L+ix) 



[L + ix) cot{kL + ikx)dx 



(22) 



decay exponentially as L ^ oo as long as we avoid the poles (i.e. kL ^ nvr). 
Choose Ltv = 4^ + ttt; then we have 



k ' 2fc' 



-a{LM+ix) 



(Lat + ix) cot{kLM + ikx)dx 



(tLm 



e '■'^'^(Ln + ix)ta.nh{kx)dx 



(23) 



This is enough to ensure exponential decay for any kL ^ nvr, since we can 
shift the right side from to — < L < Ljy + at a "cost" of less 
then jLj^e''^^'^ . Combining estimates (|2l|) and (|23| ) we see that integration 
over the part of the contour with Re 2; > decays exponentially like e"*^^^. 
In the limit L^v 00 we obtain the desired equality: 



n=l 



d 

dze-^'z—lnAiiz). 

—ioo 



(24) 



In the limit of large R, we have on the imaginary axis (Re uj = 0) 

,(3) A (4) 



Ar'Ar'iu;) 



^ sm{2ujR-ln)hi^^\uja)^[ujhi^^\uja)], (25) 



2a; 



du! 



where 



h\'\u)=ji{uj) + ini{uj). (26) 
After summation over the / index and some algebraic manipulation we obtain 



— ^ poo 

= --lmi hm^(/ + -) / duj e---c^£(/n[l- (s/e^acu))''] 



-ln[l - (sieiibiu))''] j (27) 

in agreement with (3). Note that the usual derivation of (3) involves a 
rescaling x = au. In our derivation we avoided this step and as a result, 
were able to identify the term in the cu-integral which is independent of the 
radius, and hence makes (3) divergent. 
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3 Expansion in first order 



In the previous section we reviewed the derivation of the expression for the 
Casimir energy of a conducting sphere. The Casimir energy of a dielectric 
ball, under the condition ^JeJL = \/e'fi' = 1 |]^, can be calculated at no 



further cost by repeating essentially the same steps as in the conducting case. 
Furthermore, it was shown [0, that calculations in this case, will yield the 
same result independent of the regularization used. From now on this will be 
the setting, with ^ = 1 (4) corresponding to the previous section. We now 
show that the sum over the angular index / in (3) can be carried out before 
doing the integration over the frequencies. We assume small, and expand 
the logarithm in (3) in powers of The first term in this expansion is 



oo ^ , 

-eT.(^ + 7;>jzM'- (28) 



dx 
1=1 



It turns out that the sum over / in ( j28|) can be done exactly. To this end, 
note the following identity, which can easily be obtained from the expansion 
of the Helmholtz propagator in spherical coordinates [^] using relations (5) 
and (6), 

Dix.p) = = > (2/ + l)si(xr)ei(xr)Pi(cose) (29) 

P 

^ 1=0 



p = + — 2rr cos 6 (30) 

where 6 is the angel between two vectors of lengths r and f. 

The following can, of course, be carried out for arbitrary f and r . This 
should be useful in the more general case where y^e/i 7^ v^e'/i', since the 
expressions appearing in that case involve combinations of si{x')ei{x) where 
X = y/ejiu! and x' = \J t' p'uj Jl^, but in the present case r = f = 1 is just 
what we need. Using the orthogonality relations of the P/'s 

P,{x)P,[x)dx = h,,^^ (31) 
and substituting dcos Q = —pdp we obtain the identity 

J](2/ + l)(s,eO'' = - / {^Dix,p)fpdp (32) 
1=0 •'^ 
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Thus the sum over / in (|28| ) can be carried out, as promised, and we find 



^(i - |e-4"(l - 2x)^ + 4xe-*^) (33) 

We note the appearance of a term ^ which is independent of x and thus 
causes (3) to diverge. We dispose of it by subtracting from (|33|) the appro- 
priate term corresponding to a radius b and taking the hmit b ^ oo, and 
obtain the density 

G(uj) = ^e-^'^l^l (1 + AaU + Aaltof) (34) 
Stt 

Thus, finally the first term in the expansion of (3) can be analytically 
derived. The energy to this order is 

/•oo t2 c t2 

E = / duG{u) = ^— + O(e') = ^0.0994718 + O(e'). (35) 

^ = 1 corresponds to the case of a conducting sphere. It turns out that 
this result has already been derived for a conducting sphere, by Balian and 
Duplantier |^ using a multiple scattering formalism, which may be applied 
in the case of conducting boundaries of arbitrary shape. Using a Green's 
function method, Milton, DeRaad and Schwinger evaluated the Casimir 
Energy for a conducting sphere. Numerically, to be ^0.092353. This is in 
accordance with our calculation, since the next term in the expansion of (3) 
may be evaluated using the Debye asymptotic expansion of Bessel functions 
to be —^0.007. However, for a dilute sphere (^ << 1) our result is the exact 
one. Note that our result was obtained using only subtraction of energies, 
which is in the definition of the Casimir energy. No further regularization of 
the expression (ESI) was necessary. 



4 An integral expression for higher orders in 

e 

Generalization of the method we developed in the previous section to higher 
orders in is possible. To this end, we derive a general formula for sums 
of the form Ylii^ + l/^)'^"' arbitrary power n, and coefficients a/, when 
the function f{x) = + l/2)ct/-PK^) is known. Let us first introduce the 



9 



following transform 

fix) ^ f{x,y) = - [ dtf{xy - 1/2) cost) (36) 

^ Jo 

This transform has the simple property that if f{x) = Yli <^iPi{x) then 

f{x,y) = J2^iPii^)Piiy) (37) 

I 

To see this, we write x = cos a, y = cos/9 and consider the following well 
known addition formula: 

I 

P,(cos7) = Pi{cosa)Pi{cosf3) + 2 ^ P;*^ (cos a) i^™ (cos/?) cos(m(^ - if)) 

m=l 

(38) 

where 7 is the angle between two unit vectors with spherical coordinates 

{a, 9) and (/?, v^), that is 

cos 7 = cos a cos /5 + sin a sin /3 cos(6' — ip) (39) 

setting 6' = t and ip = 0, we obtain: 

/(cos a cos /? — sin a sin P cost) = a^P^ (cos 7) = 

J2i ai [P/ (cos a) P, (cos/?) + 2 ^.^^^ P™(cosa)P^"^(cos/?) cos(mt)] (40) 

Integrating over t, we eliminate all the terms containing cost, and are left 
with '^iO,iPi{x)Pi{y) as claimed. Using (|37|) and (^) repeatedly, we obtain 

00 „i n— 1 

^(/ + l/2)af= / YldXj{f{Xi) f{xi,X2) ■ ...f{Xn-2,Xn-l)J {Xn-l)) 

[n—2)terms 

(41) 

Thus, for example, we can cast the term in the expansion of (3) in the 
form: 



327r3 



( / dx d COS Old COS 92d COS 63 / dtidt2 I I —D(x, 
^Jo J~i Jo f^\dx 



Pi) - 



(42) 

where 



pi = a/2 — 2 cos ^1 ; P2 = -\/2 — 2(cos^^i cos02 — sin^i sin^^2 costi) 

Ps = a/2 — 2 (cos 6*2 cos ^3 — sin 6*2 sin 9^ cos t2 ) ; P4 = a/2 — 2 cos 6*3 (43) 



It is not difficult to check that ( |4^) indeed converges but, unfortunately, we 
have not succeeded in integrating it analytically. 
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5 Discussion 



In this paper we developed and used a novel method appropriate for calcu- 
lations of the Casimir energy for spherical boundary conditions. Although 
the use of the asymptotic Debye approximation is very convenient, we have 
shown that it is possible to make a direct summation over angular modes us- 
ing the expansion of the Helmholtz propagator in spherical harmonics. This 
method can be applied to obtain closed integral representations for any order 
in the expansion of the Casimir energy (3). The order was explicitly 



calculated to be this result is within 0.13% of results obtained via the 



Debye approximations [jl^ , demonstrating their high accuracy. Higher orders 
in however, still remain in integral form. 

Our method should be applicable to other sums of the type ( pSf ) and 
variations on them |jlOl, which are common in the calculations for spherical 



boundaries, in order to get more regularization independent results. The fact 
that summation using angular integrals over the propagator of the Helmholtz 
equation through (^) can be done to all orders, reflects the relation between 
the mode summation technique and calculations using Green's function. This 
can be seen for the case of conducting boundaries in the multiple scattering 
formalism It would be illuminating to reveal the connections between 
our method and other methods of calculation in the dielectric case, involving 
summation over dipoles, such as calculations of the Casimir energy using 
the statistical mechanics partition function as was performed by H0ye and 



Brevik m 
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Abstract 

In this paper we sum over the spherical modes appearing in the 
expression for the Casimir energy of a conducting sphere and of a 
dielectric ball (assuming the same speed of light inside and outside), 
before doing the frequency integration. We derive closed integral ex- 
pressions that allow the calculations to be done to all orders, without 
the use of regularization procedures. The technique of mode summa- 
tion using a contour integral is critically examined. 



1 Introduction 

The Casimir effect is a remarkable consequence of QFT, exhibiting the reality 
of the zero point energy of the vacuum. This energy can be made manifest by 
studying its dependence on various constraints imposed on the electromag- 
netic field. Casimir's [1] original analysis for parallel plates yielded a force 
per unit area {h = c = 1) 



240a4 ^ ^ 

between two conducting plates separated by a distance a. The Casimir force 
for a conducting spherical shell was first calculated by Boyer [2]. The result 
was indeed puzzling: the force turned out to be repulsive. This was surprising 
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in view of the connection between Casimir's force and Van-Der-Waals forces 
which are attractive. This phenomenon may be also related to a fundamen- 
tal difference between the spherical Casimir problem and that of the plates: 
expansion of the sphere changes the distances between the points on the shell 
itself, and thus it is not clear weather the Casimir energy (associated with 
interaction among different patches on the shell) and electromagnetic the self 
energy of the shell can be separated [3]. It is not surprising, however, that 
this in turn called for a large number of subsequent calculations by various 
techniques [4, 5, 6, 7, 8]. There has been some controversy concerning the 
regularization procedures used in the calculations, bearing, for example on 
the relevance of Casimir's effect to sonoluminescence [10, 11, 12, 13]. Since 
the repulsive nature of the Casimir force for spherical boundary conditions 
is somewhat counterintuitive it seems of importance to justify it more rigor- 
ously. 

We define the Casimir energy of a conducting sphere as follows: Let 
^(a,fi) ^j^g eigenmodes of the electromagnetic field inside two concentric 
spherical shells of radii a and R {a « R), where a large sphere of radius 
R is introduced as an infrared cutoff in order to discretize the spectrum 
outside the smaller sphere as well as inside it. The zero point energy of this 
system, namely ^ iOn''^'^ is a divergent sum. The Casimir energy is then 
defined as the zero point energy difference between our configuration and a 
configuration where the inner sphere is of radius b {a « b < R) [14]^ 

Ec = lim lim lim - V (e-'^""'"''''a;i'^'^) - e"'^'"" '""'(j^'''^^) (2) 

n 

Here e^'^"" is an exponential regulator. We show that if the subtraction 
procedure is done properly, no further regularization (apart from keeping the 
naive exponential regulator, required to make the sums well defined, and 
the subtraction) is actually needed. In the sequel , we follow closely the 
procedure of mode summation using a contour integral, but our main result 
can also be applied to the Green's function method [6]. 

Most of the calculations have been carried out for the case of a conducting 
sphere and of a dielectric ball having the same speed of light inside and 
outside, that is, ^Jtjl = ji' = 1 [15], where e, are the permittivity and 
permeability of the medium surrounding a ball of radius a, and e', yu' are the 
permittivity and permeability inside the ball. It has been shown that the 

^this form of subtraction has the advantage of having a natural one to one correspon- 
dence between the eigenmodes of the constrained system (inner radius a) to those of the 
"unconstrained" one (where 6 and R are taken to oo) 



2 



Casimir energy in the latter case is given by [7, 17] 

1 poo °° 1 J 

where x = \y\,p = I + ^ and 

e = ^; s = ^ (4) 

IJi-\- IJL a 



si{x) = ixji{ix) = 




ei[x) = ixh\ [ix) = ^ Hj^yix) (6) 

and in the end of the calculation the regulator 5 is taken to zero. In the 
case ^ = 1, (3) is equivalent to the expression for the energy of a conducting 
sphere in vacuum. 

In order to evaluate (3), extensive use has been made of uniform asymp- 
totic approximations of the modified Bessel functions si{x) and ei{x). Our 
evaluation refrains from the use of these asymptotics, which, although very 
convenient in actual calculations, somewhat obscures the nature of the ex- 
pression (3). By studying the /-sum of this expression before making the 
frequency integration, using simple orthogonality considerations we arrive at 
a general formula for sums of the form ^(2/ + l)(s;e;)". We show that to first 
order in the integrand consists of a constant term (which is independent of 
the radius) and a function of a. Thus the divergence appearing in the eval- 
uation of (3) is due to integration of this constant over all frequencies. This 
divergence was already pointed out by Candelas [18], however, this constant 
term would not be present if subtraction of the Casimir energy of a large 
radius sphere is done properly. 

We start by carefully reviewing the procedure of mode summation, using 
a contour integral, as was done in some previous works [9, 7, 17]. In this 
method the final integration is carried along the imaginary axis. Justification 
for the neglect of the part of the contour away from this axis will be given, by 
concrete estimates, taking into account the existence of an infinite number 
of poles of the integrand on the real axis. 
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2 Mode summation for conducting spherical 
shells 



In order to make the discussion well defined we investigate the following 
situation [19]: A conducting sphere of radius a is placed (concentrically) 
inside a large conducting sphere of radius R. Using multipole expansion [20] 
the appropriate modes of the electromagnetic fields can be found, through 
the boundary conditions Bj. = and Ee^^ = on the shells. The eigenmodes 
are given as solutions to the equations: 
(I) Inside the inner sphere 

A\'^ = jiicoa) = (7) 
= (8) 



(21 ^ r / \i 



r—a 



II) Between the shells 

Ap^ = ji{uoR)ni{Loa) - ji{coa)ni{u;R) = (9) 



(4) 



dr 



d_ 

. dr 



rni{ior)] 



r=R 



d_ 

dr 



r=R 



dr 



[rni{ior)] 



= 
(10) 



Where equations (7), (9) refer to transverse electric modes, and (8), (10) refer 
to transverse magnetic modes. Next we define the function [19] 



AK(^;a,i?)=^^A«Ap)Af)Af) 



(11) 



A; is an analytic function, whose positive real zeros are the modes of the 
electromagnetic field inside the concentric spheres. 

In order to sum over the modes of a given I, we use the identity [7, 9] 



N 

E 



(n,l)t 



— (f dze-^'z-^lnAiiz) 



27ri 



dz 



:i2) 



Here Cn is a contour which encircles the zeros ..a;(Ar^;) of A; . Since 

there is an infinite number of modes for each /, this sum will diverge if we 
take N to infinity and a to zero. As usual, the Casimir energy should be 
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Ln + ILn 




Figure 1: A rectangular contour in the complex plane. 



identified with that obtained after the subtraction of the zero point energy 
of the unconstrained "free" system [14], 

.00 „ 

AE(a,b,R)= Urn V (/ + -) (i dze'"' zFi (13) 

N-,00 2m ^ 2 Jc^ 

Where 

Fi{z;a,b,R) = In A^^; a, i?) - ^ In A;(^; 6, i?)) • (14) 

The Casimir energy is then given by the limit: 

Ec = lim lim lim AE{a,b,R). (15) 

After the limits R ^ 00 and 6 ^ 00 are taken the remaining expression 
will have no reference to R and b. Removal of the ultra-violet cut off by 
taking the limit a ^ is done in the end of the derivation, after insuring 
convergence for all values of b and R. 

In order to make the evaluation of the contour integrals easier we take 
rectangular contours of side (Fig. 1). Our first aim is to carry the entire 
integration along the imaginary axis, by showing that the contribution of 
the part of the contour away from it vanishes as L^r goes to infinity. The 
justification of this neglect, however, is not obvious, since the integrand in 
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(13) has infinitely many singularities along the real axis, which become ever 
denser as i? ^ oo. This follows from the property that for any finite R there 
is an infinite number of field modes, and subsequently there is a corresponding 
infinite number of poles of z-j^ \n Ai{z,b, R) along the real axis. Even after 
taking the limit R ^ oo the spectrum will still consist of a discrete set 
(eigenmodes inside the inner shell) imbedded in a continuum part. Indeed, 
let us examine the behavior of A; (a;; a, R) for large arguments co. Using the 
asymptotic behavior [20] 

ji^ -sm{x-—) (16) 
n; ^ --cos(a; - y) (17) 

we obtain 

A;(a;; a, R) ^ }^ sin(2a;a - Itt) sin^(a;(i? - a)) (18) 

Thus for large arguments the modes behave as ^ or 2{R-a) '^here n is an 
integer. (Note that the density of modes for large frequencies is similar to 
that appearing in the one dimensional version of the standard Casimir effect 
for parallel conducting plates separated by a distance a inside a large box 
of length R.) Let us examine the behavior of the integrand for large values 
of \io\. Substituting the asymptotics (18), we show that the integral on a 
contour that does not pass through one of the poles, decays exponentially 
fast. The integrand is then simply: 

Fi(u; a, b, R) = 2[a{cot{2u)a — In) — cot{u){R — a))) + 
_R(cot(a.'(i? - a)) - cot(a.'(i? - b))) + 

b{cot{u {R - a)) - cot{2iob - lir))) (19) 

It is easy to check that for Im a; ^ oo and any real A, B, 

I cot{Aco) - cot{Bco) I = 0(e-2™"(^^^)i™ '^) (20) 

Thus on the upper and lower parts of the contour (Im z = ±Ljv) the integral 
is bounded by: 



.La 

-^"^^dz 



/ {iLff + z)Fi{iL]\f + z)dz < / \iL]\{ + z\12Re~ 
Jo Jo 

< 18L^i?e-^"^^ (21) 



6 



The part of the contour with Ke z = can be evaluated as follows. First 
we prove that integrals of the form 

^-a{L+ix) ^ -^^ cot{kL + ikx)dx (22) 



decay exponentially as L ^ oo as long as we avoid the poles (i.e. kL ^ mr). 



Choose L)v = 4^ + TT7-; then we have 



k '2k' 



/Ljv 
g-a(Ljv+ix)^^^ + ix) cot(A;Ljv + ikx)dx 
■Lm 



-(tLjv 



e ^'^^{Ln + ix)tanh{kx)dx 



(23) 



This is enough to ensure exponential decay for any kL ^ mr. since we can 
shift the right side from L]^ to L]^ — < L < + fj^ at a "cost" of less 
then |L|fe~°^^^. Combining estimates (21) and (23) we see that integration 
over the part of the contour with Re 2 > decays exponentially like e^"^^ . 
In the limit Lat — > cxd we obtain the desired equality: 



n=l 



(n,l)t 



nioo J 

/ d^e-'^"^— lnA;(^). (24) 

J -ioo "'^ 



In the limit of large R, we have on the imaginary axis (Re co = 0) 

A\^^A\^\u;) = ^sm{2u;R-l7r)hP{uja)^[cohP{u;a)], 
2(jj du) 



where 



(25) 



(26) 



After summation over the / index and some algebraic manipulation we obtain 
[7]: 

E, = — lim^lnn^(l + -) / duj (j£ ( /n [l - (s;e;(a(j))''] 



-ln[l - {siei{huo)Y' 



(27) 



in agreement with (3). Note that the usual derivation of (3) involves a 
rescaling x = au). In our derivation we avoided this step and as a result, 
were able to identify the term in the a;-integral which is independent of the 
radius, and hence makes (3) divergent. 



7 



3 Expansion in first order 



In the previous section we reviewed the derivation of the expression for the 
Casimir energy of a conducting sphere. The Casimir energy of a dielectric 
ball, under the condition ^Jlfi = y/FJF = 1 [7, 17] can be calculated at no 
further cost by repeating essentially the same steps as in the conducting case. 
Furthermore, it was shown [16], that calculations in this case, will yield the 
same result independent of the regularization used. From now on this will be 
the setting, with ^ = 1 (4) corresponding to the previous section. We now 
show that the sum over the angular index / in (3) can be carried out before 
doing the integration over the frequencies. We assume small, and expand 
the logarithm in (3) in powers of The first term in this expansion is 

°° 1 / 
-^'T.(^ + -2)-i^M''- (28) 

It turns out that the sum over / in (28) can be done exactly. To this end, 
note the following identity, which can easily be obtained from the expansion 
of the Helmholtz propagator in spherical coordinates [21] using relations (5) 
and (6), 

Dix,p) = = > i2l + l)siixf)eiixr)Piico^e) (29) 

p ^-^ 



p = x/r^ + — 2rr cos 6 (30) 

where 9 is the angel between two vectors of lengths r and f. 

The following can, of course, be carried out for arbitrary f and r . This 
should be useful in the more general case where ^yiJI / since the 

expressions appearing in that case involve combinations of Si{x')ei{x) where 
X = y/ejlio and x' = ^Jt' ^'uj [10], but in the present case r = f = 1 is just 
what we need. Using the orthogonality relations of the P;'s 

j'^Piix)Pjix)dx = Sij^^ (31) 

and substituting dcos 9 = —pdp we obtain the identity 

Y,{21 + l)is,e,f = - / i-Dix, p)fpdp (32) 
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Thus the sum over / in (28) can be carried out, as promised, and we find 



f (i - ie-4-(l - 2xY + 4xe-^-) (33) 

We note the appearance of a term ^ which is independent of x and thus 
causes (3) to diverge. We dispose of it by subtracting from (33) the appro- 
priate term corresponding to a radius b and taking the hmit b ^ oo, and 
obtain the density 

G{iu) = ^e-^"l'^l(l + 4a|a;| + Aa\iu\'^) (34) 

Stt 

Thus, finally the first term in the expansion of (3) can be analytically 
derived. The energy to this order is 

^ = 1 corresponds to the case of a conducting sphere. It turns out that 
this result has already been derived for a conducting sphere, by Balian and 
Duplantier [5] using a multiple scattering formalism, which may be applied 
in the case of conducting boundaries of arbitrary shape. Using a Green's 
function method, Milton, DeRaad and Schwinger [6] evaluated the Casimir 
Energy for a conducting sphere, Numerically, to be ^0.092353. This is in 
accordance with our calculation, since the next term in the expansion of (3) 
may be evaluated using the Debye asymptotic expansion of Bessel functions 
to be —2^0.007. However, for a dilute sphere (^ << 1) our result is the exact 
one. Note that our result was obtained using only subtraction of energies, 
which is in the definition of the Casimir energy. No further regularization of 
the expression (28) was necessary. 



4 An integral expression for higher orders in 

e 

Generalization of the method we developed in the previous section to higher 
orders in is possible. To this end, we derive a general formula for sums 
of the form Ylii^ + V^)^"; for arbitrary power n, and coefficients ai, when 
the function f{x) = Ylii^ + l/2)o;P;(x) is known. Let us first introduce the 
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following transform 

fix) ^ f{x,y) = - r dtf{xy - V(l-x2)(l-y2)cost) (36) 
Jo 

This transform has the simple property that if f{x) = Y2i (^iPi{^) then 

f{x,y) = J2(^iPii^)Piiy) (37) 
I 

To see this, we write x = cos a, y = cos /5 and consider the following well 
known addition formula: 

I 

P;(cos 7) = P;(cos a)P;(cos ^) + 2 ^ P;'"(cos a)P;™(cos cos{m{e - ^)) 

m=l 

(38) 

where 7 is the angle between two unit vectors with spherical coordinates 
{a, 9) and (/3, (^), that is 

cos 7 = cos a cos /3 + sin a sin (3 cos{9 — ip) (39) 

setting 6 = t and ip = 0, we obtain: 

/(cos a cos /3 — sin tt sin /3cos t) = a^P^ (cos 7) = 

Eia/[P/(cosa)PKcos^) + 2ELi^r(cos«)Pr(cos^)cos(mt)] (40) 

Integrating over t, we eliminate all the terms containing cost, and are left 
with Y^iaiPi{x)Pi{y) as claimed. Using (37) and (31) repeatedly, we obtain 

00 „i n—l 

^(l + l/2)af= / YldXj{f{Xi) f{Xi,X2) ■ ...f{Xn-2,Xn-l\ f{Xn-l)) 

(n—2jterms 

(41) 

Thus, for example, we can cast the ^'^ term in the expansion of (3) in the 
form: 

1 



327r3a 



^ / dx dcos^idcos^2C^cos^3 / dtidt2Y['i~D{x, Pi) — TT^J 

— 1 ^ -j^ 



(42) 

where 

Pi = \/2 — 2 cos 6*1 ; p2 = \/2 — 2 (cos 6*1 cos ^2 — sin 9i sin ^2 cos ti ) 

= \/2 — 2 (cos 92 cos ^3 — sin ^2 sin ^3 cos t2 ) ; p-i = \/2 — 2 cos ^3 (43) 

It is not difficult to check that (42) indeed converges but, unfortunately, we 
have not succeeded in integrating it analytically. 
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5 Discussion 



In this paper we developed and used a novel method appropriate for calcu- 
lations of the Casimir energy for spherical boundary conditions. Although 
the use of the asymptotic Debye approximation is very convenient, we have 
shown that it is possible to make a direct summation over angular modes us- 
ing the expansion of the Helmholtz propagator in spherical harmonics. This 
method can be applied to obtain closed integral representations for any order 
in the E,'^ expansion of the Casimir energy (3). The order was explicitly 
calculated to be this result is within 0.13% of results obtained via the 
Debye approximations [17], demonstrating their high accuracy. Higher orders 
in however, still remain in integral form. 

Our method should be applicable to other sums of the type (28) and 
variations on them [10], which are common in the calculations for spherical 
boundaries, in order to get more regularization independent results. The fact 
that summation using angular integrals over the propagator of the Helmholtz 
equation through (41) can be done to all orders, reflects the relation between 
the mode summation technique and calculations using Green's function. This 
can be seen for the case of conducting boundaries in the multiple scattering 
formalism [5]. It would be illuminating to reveal the connections between 
our method and other methods of calculation in the dielectric case, involving 
summation over dipoles, such as calculations of the Casimir energy using 
the statistical mechanics partition function as was performed by H0ye and 
Brevik [22]. 
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